ON THE GALILEAN INVARIANCE OF SOME DISPERSIVE WAVE 

EQUATIONS 
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Abstract. Surface water waves in ideal fluids have been typically modeled by asymptotic 
approximations of the full Euler equations. Some of these simplified models lose relevant 
properties of the full water wave problem. One of them is the Galilean symmetry, which is 
not present in important models such as the BBM equation and the Peregrine (Classical 
Boussinesq) system. In this paper we propose a mechanism to modify the above mentioned 
classical models and derive new, Galilean invariant models. We present some properties 
of the new equations, with special emphasis on the computation and interaction of their 
solitary-wave solutions. The comparison with full Euler solutions shows the relevance of 
the preservation of Galilean invariance for the description of water waves. 
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1. Introduction 

The purpose of this paper is to investigate the relevance of the Galilean invariance in 
the description of water waves. Some classical approximate models are reviewed. Those 
without the property of invariance under Galilean transformations are modified and the 
corresponding invariant versions are formulated. The new models are considered as ap- 
proximations to the full Euler equations and compared with the classical systems. The 
comparison is focused on the existence and the dynamics of solitary waves. 

A natural argument in mathematical modeling is the inheritance of the physical proper- 
ties of the phenomenon under study through the introduction of mathematical devices. In 
the case of water wave theory, the approximation to the full Euler equations leads to some 
mathematical models in which some fundamental properties of the original problem can 
be lost. This is relevant in the case of symmetries. For example, when asymptotic expan- 
sions around the still water level are performed, the invariance under vertical translations 
can be broken and the derived model is valid only in this particular frame of reference. 
Dispersive wave models possessing the property of invariance under vertical translations 
have been shown to be particularly robust for the simulation of the long wave runup, cf. 
[25]. A second symmetry, which this paper is focused on, is related to the universality of 
mechanical laws in all inertial frames of reference. The Galilean invariance (or Galilean rel- 
ativity) is one of the fundamental properties of any mathematical model arising in classical 
mechanics. This principle was empirically established by Galileo Galilei 55 years before 
the formulation of Newton's laws of mechanics in 1687, [45]. Nowadays, it is common to 
speak about this principle in terms of a symmetry of the governing equations, [47]. For 
instance, the complete water wave problem possesses naturally this property. (For a sys- 
tematic study of symmetries and conservation laws of the full water wave formulation we 
refer to Benjamin & Olver (1982), [5]). Nevertheless, numerous dispersive wave equations, 
some of them being well-known, are not invariant under the Galilean transformation. This 
issue was already addressed in a previous study by Christov (2001), [15]. We note that 
some fully nonlinear approximations such as Nonlinear Shallow Water Equations (NSWE), 
[20], improved Shallow Water Equations, [24], and the Serre or Green-Naghdi equations, 
[56, 29, 30, 37], are invariant under the vertical translation and the Galilean boost. Some 
other examples of some Boussinesq-type systems which are not Galilean invariant can be 
found in [9, 46, 6]. It is finally noted that the idea of exploiting symmetries of continuous 
equations has already been shown very beneficial in improving the behavior of underlying 
numerical discretizations, [33, 34, 14]. However, to our knowledge, the practical implica- 
tions of the loss of the Galilean symmetry are not sufficiently known. In the present study 
we try to shed some light on this issue and its influence on the approximate dispersive 
wave models. 

The paper is organized as follows. In Section 2 we review some classical models in water 
wave theory and present the invariant counterparts of the BBM equation and the classical 
Peregrine system. A way to asses the accuracy of the invariantized models is by comparing 
their travelling wave solutions to the corresponding ones of the full Euler equations and of 
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the rest of the classical models under consideration. In Section 3, solitary- wave profiles of 
the new models are generated numerically. The comparison with the wave profiles of other 
approximations (some with exact formulas) and of the Euler equations (by using Tanaka's 
algorithm and Fenton's asymptotic solution) is established in terms of the amplitude-speed 
and amplitude-shape relations. The interactions of solitary waves for the new models, with 
head-on and overtaking collisions, are studied in Section 4. Finally, the main conclusions 
of this study are outlined in Section 5. 

2. Mathematical models 

Consider an ideal fluid of constant density along with a cartesian coordinate system 
in two space dimensions (x,y). The y-axis is taken vertically upwards and the a:-axis 
is horizontal and coincides traditionally with the still water level. The fluid is bounded 
below by an impermeable horizontal bottom at y — —d and above by an impermeable free 
surface at y = r](x, t). We assume that the total depth h(x, t) = d + r](x, t) remains positive 
h(x, t) ^ ho > at all times t. 

Assuming that the flow is incompressible and irrotational, the governing equations of 
the classical water wave problem, [58], are the following: 

<Pl + <f>yy = -d ^y^ V {x,t), (2.1) 

rjt + 4>xVx - 4> y = y = r)(x,t), (2.2) 

<f>t + \(<Pxf + \{<t>y? + 9V = y = v(x,t), (2.3) 

0, = y = -d, (2.4) 

with being the velocity potential (by definition, the irrotational velocity field u = <p x and 
g the acceleration due to the gravity force. The water wave problem possesses Hamiltonian, 
[53, 61, 11, 55], and Lagrangian, [43, 16], variational structures. 

The symmetry group of the complete water wave problem (2.1) - (2.4) was described 
by Benjamin & Olver (1982) in [5]. In particular, the full formulation of the water wave 
equations admits the Galilean boost symmetry and the invariance under the vertical trans- 
lations (the latter issue will be addressed by the authors in a future work). However, the 
water wave theory has been developed from the beginning by constructing various approxi- 
mate models which may conserve or break some of the symmetries, [19]. Below we consider 
several classical models and discuss their Galilean invariance property. 

2.1. The KdV equation. The unidirectional propagation of long waves in the so-called 
Boussinesq regime (where the nonlinearity and the dispersion are of the same order of 
magnitude), [10, 6], can be described by the celebrated Korteweg-de Vries (KdV) equation 
[35, 32], which in dimensional variables is written in the form: 

/—, 3 d2 r~r 

u t + \/gdu x + -uu x + —\/gdu xxx = 0, (2.5) 

2 6 

where u(x,t) is the horizontal velocity variable which is usually defined as the depth- 
averaged velocity, [52], or the fluid velocity measured at some specific water depth, [9]. 
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Some very well-known properties of (2.5) are reviewed (see e.g. [28, 38, 44]). First, the 
KdV equation is an integrable model with the following two-parameter family of solitary 
wave solutions: 

u(x, t) = uosech 2 (^k(x — c s t — xo)) , c s = \fgd-\ — -, (nd) 2 = ° , m > 0, xq G K. 

2 y/gd 



The initial value problem of the Korteweg-de Vries (KdV) equation possesses a Hamil- 
tonian structure 

6H 

ou 

(where 5 denotes the variational derivative) in a suitable phase space of functions vanishing, 
along with some of their derivatives, at infinity. The skew-symmetric operator J and the 
Hamiltonian functional H are 

d 2 



—d x , H = — I \fgdu 2 + |w 3 — \j~g~d 
2 Jn L 



dx. 



The Hamiltonian H is the third conserved quantity of the well-known hierarchy of invariants 
for (2.5), [38]. 

The central question in our study is the Galilean invariance of model equations. We can 
show that the KdV equation (2.5) possesses this property. The procedure is as follows. 
We choose another frame of reference which moves uniformly, for example, rightwards 
with constant celerity c. This symmetry is expressed by the following transformation of 
variables: 

x x — |ct, t — > t, u — » u + c. (2.6) 
In this moving frame of reference (2.5) becomes: 

3 f— 3. , d 2 r- 

Ut - t; cu x + V9du x + -{u + c)u x + —\Jgdu xxx = 0. 
2 2 o 

After some simplifications one can recover the original KdV equation, which completes the 
proof of the invariance. 

In order to assess the relative magnitude of various terms in equation (2.5) scaled vari- 
ables are introduced. The classical long wave scaling is the following: 

x':=- y':= y - t' := U, rj' := A u' := -=, (2.7) 
t d d a y/gd 

where h , a, £ are the characteristic water depth, wave amplitude and wave length re- 
spectively. Using these three characteristic lengths we can form three following important 
dimensionless numbers: 

*:= 5 , ^:=( ? ), 3 = ^. (2.8) 

Parameters £ < 1 and fi 2 1 characterize the wave nonlinearity and dispersion, while 
the so-called Stokes number S measures the analogy between these two effects. In the 
Boussinesq regime the Stokes number is of order of one, S = 0(1), which establishes that 
the dispersion and the nonlinear effects are comparable. 
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Using these dimensionless and scaled variables the KdV equation (2.5) can be written 
in the form: 

3 ^ 
u t + u x + -euu x + —u xxx = 0. 

2 o 

where the primes have been dropped. Formulas (2.7) and (2.8) will also be used in some 
of the developments below. 

2.2. The BBM equation. Benjamin, Bona & Mahony (1970), [4], (see also [51]) proposed 
the following modification of the KdV equation, known as the BBM equation: 

r~r 3 d 2 , . 

u t + \/gdu x + -uu x - —u xxt = 0. (2.9) 
z o 

The main idea for the derivation of this model is to use the lower order relation between 
time and space derivatives in order to modify the higher-order dispersive termut = — u x + 
0(e + v 2 ) which implies u xxx = -u xxt + 0(e + /i 2 ). 

One of the main practical motivations for this modification is to improve the dispersion 
relation properties of the KdV equation. Specifically, unlike the KdV equation, the phase 
and group velocities of the BBM equation have a lower bound. It is also noted that the 
BBM equation has the following solitary wave solutions: 

u(x, t) = u sech 2 (i/t(a; — c s t — x )), c s —y/gd+—, (nd) 2 = 1 — , x e R. 

2 y gd+ 2% 

The BBM equation is not integrable, but it can also be written as an infinite-dimensional 
Hamiltonian system 

5U 

u t = J—, 
ou 

where the operator J and the Hamiltonian functional H, are defined as: 

J = (1 - %d xx Y x ■ (-d x ), U = [\fgdu 2 + lu 3 ] dx. (2.10) 

and the structure is defined on a phase space similar to that of the KdV equation. 

As far as the Galilean invariance is concerned, the change of variables (2.6) applied to 
(2.9) leads to 

3 r- 3, , d 2 d 2 

u t - -cu x + \/gdu x + -{u + c)u x - —u xxt + —cu xxx = 0, 
2 2 o 4 

and after some algebraic simplifications we obtain: 

r- 3 d 2 d 2 
u t + \/gdu x + —uu x - —u xxt + —cu xxx = 0. 
2 o 4 

Since there is at least one new term (^cu xxx ) appeared in the previous moving frame of 
reference, the BBM equation is not Galilean invariant. The relevance of this drawback 
always puzzled the researchers, cf. [15]. 
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2.3. The iBBM equation. We now propose a modification to the classical BBM equation 
which allows us to recover the Galilean invariance property. Furthermore, the idea behind 
the arguments below can be extrapolated to other models. The strategy is to add a new 
term which will vanish the non-invariant contribution of the BBM dispersion u xx t under 
the transformation (2.6). The resulting equation, which will be called invariant Benjamin- 
Bona-Mahony (iBBM) equation, takes the form: 

/— 3 d 2 d 2 . N 

u t + ygdu x + -uu x - —u xxt - —uu xxx = 0. (2.11) 
2 6 4 

It is straightforward to see that (2.11) is invariant under the Galilean transformation (2.6). 

The modification proposed above becomes more transparent in scaled variables. The 
application of the long wave limit (2.7) to (2.11) leads to 

3 /x 2 e/i 2 

Ut + u x + -euu x - —u xxt —uu xxx = 0. 

2 6 4 

One can observe that the last term on the left hand side, responsible of the Galilean invari- 
ance of (2.11), is a nonlinear term of order 0(en 2 ) and consequently, it is asymptotically 
negligible in the BBM formulation. Since this additional term is nonlinear, the linear 
dispersion relation of (2.9) is not modified. Its effect will be studied thoroughly in the 
following sections. 

Remark 1. Unlike the BBM equation (2.9), this invariant version (2.11) does not possess, 
to our knowledge, a Hamiltonian structure. However, it is possible to propose an invari- 
antization which preserves this variational formulation along with the Galilean invariance. 
The alternative given by the equation 

u t + \fgdu x + \uu x - ^—u xxt - (2u x u xx + uu xxx ) = 0. (2.12) 
2 6 4 

has an additional higher-order nonlinear term which allows for a non-canonical Hamilton- 
ian structure. In this case, the operator J = (1 — ^d^)^ 1 ■ {—d x ) is the same as for the 
BBM equation and the Hamiltonian H is 



« = i 



E 



1 



2 1 q d 
gdu + -u + —uu r 
2 4 



dx. 



We underline some similarity between equation (2.12) and several proposed earlier models 
such as the Camassa-Holm [12], Burgers-Poisson [27] and Degasperis-Procesi [21] equa- 
tions. 

We now look for travelling wave solutions of (2.11) of the form: 

u(x, t) = w(£), £ := x — c s t, (2-13) 

where c s is the solitary wave speed. We also assume that w(£) decays to zero along with 
all derivatives when |£| — > oo. Substituting (2.13) into the iBBM equation (2.11) leads to 
the ordinary differential equation: 

3 d 2 d 2 

gd - c s )u' + -(u 2 )' + cs— u'" - —uu'" = 0, (2.14) 
4 6 4 
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Figure 1 . Phase plane of the invariant models with c s > y/gd: (a) (iBBM) . 
(b) (iPer). In both cases, the trajectory — > A — > B — )■ represents a 
solitary wave. 



where the prime denotes differentiation with respect to £. Using the boundary conditions 
at the infinity, the identity uu" 1 = (uu" — j(u') 2 )', and after one integration, equation (2.14) 
becomes: 

J2 j2 



(y^d - c> + |u 2 + c s ju" - j (±(u') 2 - uu") = 0, (2.15) 
that can be written ClS db system 

u' = v, (2.16) 

^r|y(( c «-^"-i" 2 -¥" 2 )' (m) 



v' 



Now it can be checked that when c s > ygd, the origin u = v = is a saddle point, 
as depicted in Figure 1(a), which shows the corresponding phase plane. The homoclinic 
trajectory O — > A — > B — > O represents a solitary wave. (The MATLAB code for this 
figure can be found in [48].). 

2.4. The Peregrine system. Under the assumptions described above, D.H. Peregrine in 
1967 derived the following system of equations governing the two-way propagation of long 
waves of small amplitude in the Boussinesq regime [52]: 

Vt +((d + r ] )u) x = 0, (2.18) 

d 2 , s 

u t + uu x + gr) x - —u xxt = 0, (2.19) 

where u(x,t) is now defined as the depth averaged fluid velocity, r](x,t) is the deviation 
of the free surface of the water from its rest position. This system is also known as the 
Classical Boussinesq system, cf. [6] and will be denoted by (cPer). In [49], the existence and 
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some properties of solitary wave solutions of (2. 18)- (2. 19) are obtained, without explicit 
formulas. On the other hand, to our knowledge, a Hamiltonian structure has not been 
found, [6]. 

We now study the Galilean invariance of (2.18)-(2.19). In this case, the Galilean trans- 
formation takes the following form: 

x — > x — ct, t — > t, T] — )■ Tj, u — > u + c. (2.20) 

The mass conservation equation (2.18) in new variables reads: 

Vt~cr] x + ((d + r})(u + c)) x = 0. 

After simplifications one can see that this equation remains invariant under the transfor- 
mation (2.20). 

Now let us consider the momentum balance equation (2.19). In the moving frame of 
reference this equation becomes: 

d 2 d 2 

Ut - cu x + (u + c)u x + gr] x - —u xxt + c—u xxx = 0, 
and after some manipulations 

d 2 d 2 
u t + uu x + gr] x - —u xxt + c—u xxx = 0. 

As in the BBM case, a new dispersive term showing that the system 

(2.18) -(2.19) is not Galilean invariant. 

2.5. The iPeregrine system. Following the same technique as in the case of the BBM 
equation, it is possible to propose a modification to the classical Peregrine system (2.18)- 

(2.19) which will allow us to recover the Galilean invariance property. It can be done in a 
similar way leading to the iBBM equation (2.11). The corresponding system reads: 

Vt +((d + V )u) x = 0, (2.21) 

d 2 d 2 . , 

u t + uu x + gi] x - —u xxt - —uu xxx = 0. (2.22) 

Note that since the mass conservation equation is invariant, it is not modified in the new 
version. Now it is straightforward to check the invariance of equation (2.22). Therefore, 
system (2.21)-(2.22), which will be called invariant Peregrine system or (iPer) for the sake 
of conciseness, is Galilean invariant. In dimensionless and scaled variables, the system 
reads: 

r] t +({l+er})u) x = 0, 

fi 2 efi 2 
u t + euu x +rj x - —u xxt —uu xxx = 0, 

and one can see that the new term is of higher-order and, asymptotically speaking, neg- 
ligible. As in the case of the Peregrine system, equations (2.21)-(2.22) do not possess, to 
our knowledge, a Hamiltonian structure. 
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Finally, we can look for travelling wave solutions of system (2.21)-(2.22) of the form 

T](x, t) = 77(0, u(x, t) = £ := x - c a t, 

where 77 and u decay to zero, along with their derivatives, as |£| — > 00. After substituting 
this representation into the governing equations (2.21)-(2.22) they become: 

-c s T]'+ ((d + r])u)' = 0, 

1 d 2 d 2 

-c s m + -{u ) + grj +c s —u - —uu = 0. 

An integration of the mass conservation equation and the decay at infinity lead to 

u = 4^~- (2.23) 
d + r] K ' 

Then the momentum balance equation can be integrated once and we have 

Similarly to the case of the (iBBM), one can see that, (2.24) written as a first order system 
and when c 2 — gd> 0, the origin is a saddle point; the phase plane sketched in Figure 1(b) 
also shows a solitary wave, in the form of a trajectory O — > A — > B — > O. 

2.6. The Serre equations. In order to complete the presentation of our model equations, 
we consider the fully-nonlinear system referred to as the Serre, [56, 3]: 

h t + [hu] x = 0, (2.25) 
u t + uu x + gh x = I hT x [h 3 (u xt + uu xx - u 2 x )] x . (2.26) 

Solitary wave solutions of (2.25)-(2.26) are explicitly known. They are given by the for- 
mulas: 

7](x, t) = aosech 2 (\k(x — c s t — x )) , u = ~r~ — > (2.27) 

d + rj 



where c s = \/ g(d + Oq) and (nd) 2 = J^-- 

As pointed out by Li, [39] and [40], the Serre equations possess a Hamiltonian structure: 

jit) \5n/5q / 
where the Hamiltonian functional H and the operator J are given by 

o K 

hd x q x + qd x J ' 

The variable q is sometimes referred to as the potential vorticity flux and is defined by 

q ■= hu - \ [h 3 u x ] x . 

The Serre equations (2.25)-(2.26) can be shown to have the Galilean invariance property. 
For the mass conservation equation (2.25) we refer to Section 2.5. Thus it remains to check 



H — g 



E 



hu z + \h 3 u 2 x + grj 



dx, 
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this property for the momentum conservation equation (2.26). If we make the change of 
variables t — )■ t, x — >■ x — ct, h — )■ h and u — > u + c as before, equation (2.26) becomes: 

u t - cu x + (u + c)u x + gh x = \ h~ l [h 3 (u xt - cu xx + (u + c)u xx - u 2 x )] x , 

and after two simple algebraic simplifications one can recover the original equation (2.26). 

3. Numerical computation of travelling waves 

In the previous section we presented several classical models arising in water wave theory. 
Moreover, we proposed two novel equations, the iBBM equation and the (iPer) system, 
with the aim of incorporating the property of invariance under the Galilean transformation, 
lost by the original BBM equation and the (cPer) model. The purpose of this and the next 
sections is to compare these models through the computation of their respective solitary 
wave solutions and, whenever possible, with the solitary waves of full Euler equations. 
Some of the models, such as the KdV, BBM and Serre equations, possess explicit formulas 
for these solutions, while the iBBM equation, the (iPer) system and the (cPer) system 
do not. Consequently, the latter have to be constructed numerically. In the next section 
we give the great lines of the numerical procedure to this end, as well as the method we 
use to construct approximate solitary waves of the Euler equations. Both techniques will 
be applied in the comparative study. We note that in all cases we will use the models in 
nondimensional but unsealed variables (taking e = /i 2 = 1). 

3.1. Computation of travelling-wave profiles. The Petviashvili method. The in- 
vestigation for travelling wave solutions in one-dimensional systems typically leads to a set 
of differential equations of the form 

CU = N(U), (3.1) 

for some differential operators C (linear) and M (nonlinear). The numerical resolution of 
the preceding system can be done in many different ways (see [59] and the references therein 
as a modest representation of the related literature). Among all the possibilities, the so- 
called Petviashvili method will be used in our computations. This method stems from 
the pioneering work of V.I. Petviashvili (1976), [54]. It is based on a modification of the 
classical fixed point iteration (which in these cases is usually divergent) and it is formulated 
as follows. Given an initial profile Uq, the Petviashvili iteration generates approximations 
U n of the original solution of (3.1) following the formulas 

CU n+1 = M^N(U n ), (3.3) 

where (•, •) denotes the usual L 2 — norm and 7 is a free parameter that controls the con- 
vergence of the method. The term (3.2) is called the stabilizing factor. See [50, 36] for 
details, generalizations and some local convergence results. 

In this study, (3.2)-(3.3) is applied to compute solitary wave profiles in the following 
cases and with the corresponding operators, namely: 
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(a) invariant Benjamin-Bona-Mahony (b) invariant Peregrine (iPer) system 
(iBBM) equation 

Figure 2. Solitary wave profiles for different speeds c s . 

• (iBBM): Cu = (VP - c,)u + c,£u", Af(u) = f - uu") - f « 2 . 

• (cPer): £ M = c s (u - f «") , = £ + 

. (iPer): £ M = c a (u - f u") , JV(u) = T + ^- y(^- uu ") > 
where c s is the solitary wave speed. 

3.1.1. Generation and numerical evolution of the profiles. In many cases, the method (3.2)- 
(3.3) can be efficiently implemented by using Fourier techniques, [50, 36]. Specifically, our 
implementation for the three systems has been performed by considering the corresponding 
periodic problem and using a pseudospectral representation for the approximations to the 
profiles. As an initial iteration, a solitary wave solution (2.27) of the Serre equations or 
the third-order asymptotic solution of Grimshaw, [31], can be considered. The iterative 
procedure is continued until the difference between two consecutive iterations in the 
norm, or the norm of the residual is less than a prescribed small tolerance which, in 
our case, is of order O(10~ 15 ). The convergence is reached in 10-20 iterations. In order 
to illustrate better the transformations that solitary waves undergo while we gradually 
increase the propagation speed parameter c s , we superpose several profiles on the same 
Figure (Figure 2(a) corresponds to invariant Benjamin-Bona-Mahony (iBBM) system and 
Figure 2(b) to invariant Peregrine (iPer)) system. 

In order to assess the accuracy of the computations, the three models have been nu- 
merically integrated in time using the computed solitary wave profiles as the initial con- 
ditions. Some error indicators measuring the accuracy of the numerical approximation 
of the solitary waves have been computed. The numerical method for the corresponding 
initial-periodic boundary value problem consists of a pseudospectral- Galerkin method for 
the semi-discretization in space and the classical, explicit fourth-order Runge-Kutta scheme 
for the time integration, [23]. 




Figure 3. Some error indicators for the propagation of a solitary for the 
Peregrine and the invariant Peregrine systems 

In the experiments below, we study the propagation of a solitary-wave profile, generated 
by (3.2)-(3.3) with speed c s = 1.1 in the interval [-128, 128], with N = 2048 points for the 
pseudospectral approximation and and spatial and time step sizes dx = 1.25 x 10~ 2 ,<it = 
1.25 x 10 -3 respectively. The results correspond to (cPer) and (iPer) systems. The (iBBM) 
case has also been implemented, with similar results. 

We study the evolution of two parameters: the normalized amplitude error and the shape 
error. The first one is computed by comparing, at each time level, the initial amplitude 
of the profile (generated by the Petviashvili method) with the corresponding amplitude of 
the numerical solution. Their computation is implemented as in e. g. [22]. Figure 3(a) 
shows the temporal evolution of this amplitude error up to a final time T = 100, for both 
(cPer) and (iPer). We observe that for the specific values of dx and dt the amplitude is 
conserved up to 10 decimal digits in both cases. 

On the other hand, the shape error is computed by comparing the numerical solution 
with time translations of the initial profile with the prescribed speed and minimizing the 
differences (see [23] for the details). The results displayed in Figure 3(b) show a virtually 
constant evolution of this error, which in both cases is of order O(10~ 10 ). These results 
confirm the accuracy of the technique used to generate the solitary wave profiles and of the 
numerical code for the time evolution. The latter will be used for the experiments below. 

3.2. Solitary waves of the Euler equations. The construction of approximations to 
travelling wave solutions for the 2D Euler equations with free surface (to be considered 
as reference solutions for the approximate models) will be performed with two techniques. 
One by using the Tanaka's algorithm, [57], and the other by using asymptotic solutions 
given in, [26, 42]. 

3.3. Numerical results. In this section we compare the solitary waves of the following 
models: 
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• The Korteweg-de Vries (KdV) equation (2.5). 

• The Benjamin-Bona-Mahony (BBM) equation (2.9). 

• The invariant Benjamin-Bona-Mahony (iBBM) equation (2.11). 

• The Peregrine classical Peregrine (cPer) system (2.18), (2.19). 

• The invariant Peregrine invariant Peregrine (iPer) system (2.21), (2.22). 

• The Serre equations (2.25), (2.26), 

from either the analytical formula (when possible) or the computations with the Petvi- 
ashvili method. The comparison is established between them and with those of the full 
Euler equations with free surface. The latter are computed by two ways: from the Fenton 
asymptotic ninth order solution and the Tanaka's solution of the two-dimensional incom- 
pressible Euler equations with free surface. The study is focused on the amplitudes and 
shapes of the computed free surface elevation ??(£), provided by the models for the same 
prescribed value of the propagation speed parameter c s . 

3.3.1. Solitary wave speed-amplitude relation. Figure 4(a) shows an amplitude-wave speed 
diagram for the models considered. First, an approximate relation between the amplitude 
and speed for the solitary waves of the full Euler system is computed, by using the Tanaka's 
algorithm, [57], and Fenton's 14th order, [42]. They virtually give the same results and 
these are compared with the relation obtained by each of the models. We observe that the 
Serre equations (and, consequently, the other systems, that have weaker nonlinearities) are 
known to provide a relatively good approximation to solitary wave solutions of full Euler 
equations in a range of amplitudes not greater than 0.5, cf. e.g. [41, 13]. Therefore, our 
attention is focused on solitary waves with these amplitudes, as it is observed in Figure 
4(a). Figure 4(b) shows a magnification for the largest amplitudes. 

Corresponding to the comparison between the solitary waves of the models at hand, we 
note that the non Galilean invariant models, the Benjamin-Bona-Mahony (BBM) equa- 
tion and the classical Peregrine (cPer) system, tend to underestimate the solution speed 
for a given amplitude. On the other hand, the curves corresponding to their invariant 
counterparts along with the fully nonlinear Serre equations lie above the reference solu- 
tion. In particular, the results for the invariant Benjamin-Bona-Mahony (iBBM) equation 
are very close to those of Serre equations. A surprising fact is that the amplitude-speed 
relation given by the invariant Peregrine (iPer) system is superposed with that of the Serre 
equations up to the graphical resolution. 

3.3.2. Solitary wave shapes. A second comparison between the shape of the computed 
solitary waves, is presented in Figures 5(a)-5(f). They illustrate the cases of small (a/d ~ 
0.1, see Figure 5(a) and a magnification on the wave crest in Figure 5(b)), moderate 
(a/d « 0.22, see Figures 5(c) and a magnification in Figure 5(d)) and large (a/d m 0.4, see 
Figures 5(e) and 5(f)) solitary wave amplitudes of the models (within the range mentioned 
above). According to these results, it is observed that the invariant Benjamin-Bona- 
Mahony (iBBM) equation and the invariant Peregrine (iPer) system approximate much 
better the amplitude of the reference solution (represented, in this case, by the Tanaka's 
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Figure 4. The curves corresponding to the iPeregrine system and the Serre 
system are superposed up to the graphical resolution. 



solution) than the non-Galilean invariant counterparts, and they stay very close to the 
results of the Serre system near the crest. 

As a measure of the level of approximation to solitary wave solutions of the Euler system, 
these and the previous results show the benefits of taking into account the invariantization 
process in the approximate models. 



4. Solitary waves' dynamics 

We complete the numerical experiments by studying the effects of the Galilean invariance 
property in the evolution of solitary waves. Specifically, we first compare, by numerical 
means, head-on collisions of two solitary waves of the invariant BBM and Peregrine equa- 
tions with those of their corresponding not invariant models. Then we will study the 
evolution of the four systems when when a solitary wave profile of the Euler equations is 
taken as initial condition. 

It is known that the tails produced by the interaction of two solitary waves are sensitive to 
both the linear terms (characterizing the linear dispersion relation) and the nonlinearities. 
For example, two solitary waves of the KdV equation interact in an elastic way without 
producing dispersive tails at all, [60], while the collision of solitary waves of the BBM 
equation will produce dispersive tails and probably small- amplitude nonlinear pulses as 
the main indication of an inelastic interaction, [7, 8]. In the new Galilean invariant models, 
the new nonlinear terms are of order Efi 2 and their effects on the interaction will be studied 
here. The same code introduced in Section 3 is used for the numerical computations, as well 
as the Petviashvili method (3.2)- (3.3) to generate solitary wave profiles when necessary. 
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Figure 5. Solitary waves of different speeds c s 



16 A. DURAN, D. DUTYKH, AND D. MITSOTAKIS 

1.014 
1.012 
1.01 

1.008 

K 

1.006 
1.004 
1.002 



1.1 1.15 1.2 1.25 1.3 1.35 1.4 

Speed c s 

Figure 6. Symmetric head-on collision. Ratio K of amplitudes before and 
after the collision, for Peregrine and invariant Peregrine systems. 



(a) Peregrine system (b) iPeregrine system 



C s 




A a fter 




c s 






1.1 


0.2177418 


0.2177417 




1.1 


0.21 


0.209978 


1.2 


0.4757297 


0.4757202 




1.2 


0.44 


0.439365 


1.3 


0.7822906 


0.7821674 




1.3 


0.69 


0.686027 


1.4 


1.1476304 


1.1469096 




1.4 


0.95 


0.946593 



Table 1. Symmetric head-on collision. 



4.1. Head-on collisions of solitary waves. A first group of experiments concerns head- 
on collisions. The classical Peregrine (cPer) and the invariant Peregrine (iPer) systems 
have been considered, by constructing, in both cases, two solitary waves on the interval 
[—256,256] with speeds c 8> \ = 1.15 and c Sj 2 = 1-05 (translated appropriately such as their 
amplitudes achieved on —50 and 50 respectively) and travelling in opposite directions. 
These solitary waves are of small amplitude and their shapes are almost the same for both 
models. The code uses N = 4096 nodes, a spatial step size of 1.25 x 10 _1 and 5 x 10~ 3 as 
the time step. 

Figure 7 shows the t]— evolution of the head-on collision for both models and at several 
times. A tail behind each tallest wave after the collision is observed. This is larger in the 
case of the invariant Peregrine system (of the order of 10~ 4 ) than in the case of the classical 
Peregrine system (approx. 10~ 4 ), see Figure 7 (d) and (f). During the collisions, a similar 
phase shift also takes place, see Figure 8. 

A final comparison is established in terms of the degree of inelasticity of the interaction. 
This can be measured by using several parameters, [1, 17]. In each symmetric 
head-on collision has been implemented; that is, two solitary waves with the same spped 
c s travelling in opposite directions. After the interaction, both solitary waves emerge with 
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Figure 7. Head on collision of two solitary waves for the Peregrine and 
invariant Peregrine systems. 



similar amplitudes A a f ter , but below the initial one A^u, as is observed in Table 1. Then 
a ratio K of the amplitude of the waves after the collision to their amplitude before the 
collision has been computed. Figure 6 shows the behavior of this value, as a function of the 
speed parameter c s and for both systems. The results reveal a higher inelastic collision in 
the case of the invariant system, in accordance to what is observed in the full water wave 
model [18]. 
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(a) (b) 



Figure 8. Phase diagram of the overtaking collision: (a) Peregrine system, 
(b) Invariant Peregrine system. 

4.2. Overtaking interactions of solitary waves. In order to study the overtaking colli- 
sion of solitary waves, we first consider the classical and the invariant Peregrine systems and 
construct, for both, two solitary waves on the interval [—1024, 1024] with speeds c S]1 = 1.15 
and c s 2 = 1-05 (translated appropriately such as their amplitude achieved on —50 and 50 
respectively) travelling in the same direction, with N = 16384 and time step of 5 x 10~ 3 . 
In Figure 9 we observe that the basic characteristics of the interaction are similar for both 
systems, (cf. also [2]). Specifically, the interaction is again inelastic; after the collision a 
tail, of apparent dispersive nature, behind the smallest wave (moving to the right) and a 
small N-shape wavelet (moving to the left) are generated. Moreover, a small phase shift 
and a change in shape can be observed in the solitary pulses (Figure 8). The wavelet 
generated in the case of the classical Peregrine system has an inverse N-shape while the 
invariant version has an N-shape. 

With the same input data, an overtaking collision for the case of the BBM and iBBM 
equations is shown in Figure 10 at several times. In this case, on the contrary, no N-shape 
wavelet is observed and only a tail behind the waves appears. 

We can conclude that these two groups of results show that, as expected, the (9(e/x 2 ) 
nonlinear terms of the invariant models do not change significantly, in a qualitative sense, 
the behaviour of the interactions of solitary waves provided by the corresponding not 
Galilean invariant system. 

4.3. Comparison with Euler equations. Finally, we study the evolution of a solitary 
wave of the Euler equations when we use it as initial condition to the approximate models. 
Specifically, we consider an approximate solitary wave $h(x) of amplitude A = 0.2 (and 
speed c s ~ 1.095490471188718) obtained by using Fenton's ninth order asymptotic solution. 
In the case of the Peregrine and the iPeregrine systems we use for initial velocity Uq(x) = 
c s rj (x) / (d + r] (x)) . We remind that due to the mass conservation property this formula is 
exact (see equation (2.23) for the velocity given the surface elevation r\ for both Peregrine 
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Figure 9. Overtaking collision of two solitary waves for classical Peregrine 
(cPer) and invariant Peregrine (iPer) systems. 



and iPeregrine systems). In Figure 11 we present the solution at T = 100. We observe 
that the initial condition is resolved into a new solitary wave followed by a dispersive tail. 
In the case of the classical models, the dispersive tails appear to be smaller. Figure 12 (a) 
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Figure 10. Overtaking collision of two solitary waves for Benjamin-Bona- 
Mahony (BBM) and invariant Benjamin-Bona-Mahony (iBBM) models. 



shows the shape error of the solution (i.e. how much different is the solution for being the 
exact solitary wave of the Euler equations) while Figure 12 (b) presents the amplitude of 
the solution as a function of time t. From these two figures we observe that in the case 
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FIGURE 11. Evolution of a solitary wave of the Euler equations, with 
(BBM), (iBBM), (cPer) and (iPer) models. (A = 0.2) 



of classical models the emerging solitary waves are closer in shape and amplitude to the 
original solitary wave solution of the Euler equations than the respective solitary waves of 
the invariant models. 



5. Conclusions 

In the present work the influence of Galilean invariance in several equations arising in 
water wave modelling is studied. We propose the modification of some existing not invari- 
ant models in order to include this fundamental property. The technique introduced here 
consists of adding higher-order terms from the approximation of the governing equations. 
These terms are asymptotically negligible and consequently, the modified models are still 
valid in the appropriate regime. As a case study, corresponding modifications of two not 
invariant models, the Benjamin-Bona-Mahony (BBM) equation and the classical Peregrine 
system, are presented. The comparison with reference solutions to the full Euler equations 
shows that this extra-term is beneficial for the description of the travelling wave solutions 
in several ways. First, this modification improves the SW amplitude-speed relation which 
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Figure 12. Shape error and amplitude of the solution as a function of time 



lies closer to the Tanaka's and Fenton's solutions. In this regard, we obtain a surpris- 
ing performance of the iPeregrine system (2.21), (2.22) with the amplitude-speed relation 
undistinguishable from the fully- nonlinear Serre equations (2.25), (2.26). Moreover, the 
amplitudes of solitary wave solutions to the invariant models are closer to the correspond- 
ing full Euler solutions than classical counterparts. The comparison is finished off with a 
numerical study of head-on and overtaking collisions. Compared to the behaviour observed 
in the not invariant equations, the higher order nonlinear terms incorporated in the new 
models do not affect qualitatively the inelastic character of the interactions. However, a 
relevant difference in the degree of inelasticity is observed, being higher in the case of the 
invariant versions. This behaviour is closer to what has been observed in the case of the 
full Euler equations. 
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